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Abstract 

• Conventional nucleon-nucleon potentials with strong short-range repulsion require 

contributions from high-momentum wave function components even for low-energy 
observables such as the deuteron binding energy. This can lead to the misconcep- 
' tion that reproducing high-energy phase shifts is important for such observables. 

Interactions derived via the similarity renormalization group decouple high-energy 
and low-energy physics while preserving the phase shifts from the starting potential. 
They are used to show that high-momentum components (and high-energy phase 
shifts) can be set to zero when using low-momentum interactions, without losing 
. information relevant for low-energy observables. 
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1 Introduction 



High-momentum degrees of freedom do not automatically decouple from low- 
energy observables, especially for conventional nucleon-nucleon (NN) poten- 
tials. For instance, nuclear forces are typically fit to scattering data up to where 
inelasticities start to become significant, ~ 350 MeV or relative momenta 
k ~ 2fm~ . However, most NN potentials have significant high-momentum 
(k > 2fm _1 ) off-diagonal matrix elements that require summations over high- 
energy intermediate states, even if one is calculating low-energy observables. If 
such interactions are simply truncated at 2fm~ , the deuteron binding energy 
along with S-wave phase shifts down to zero energy are drastically altered. 
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This can lead to the misconception that details of strong-interaction dynam- 
ics above some energy scale are relevant to low-energy nuclear structure and 
reactions. Such a brute-force cutoff, however, does not disentangle high-energy 
(short- distance) features from low-energy (long-distance) observables. To do 
so, it is necessary to integrate out (and thus separate) irrelevant short- distance 
details from their effects on low-energy observables. This is achieved by the 
renormalization group. 

Renormalization group (RG) transformations that lower a cutoff in relative 
momentum have been used to derive NN potentials that have vanishing matrix 
elements for momenta above the cutoff. Such interactions, known generically 
as Viowfc, show greatly enhanced convergence properties in nuclear few- and 
many-body systems for cutoffs of order A = 2fm~ 1 or lower [Tll2T3llll5ll6lr7] . 
However, the initial NN potentials (typically cut off at 4-5fm _1 ) predict non- 
zero elastic phase shifts to much higher energies than ~ 350 MeV, which in 
some channels are semi-quantitatively consistent with experiment. In contrast, 
phase shifts from Vi ow fc as usually implemented are zero above the cutoff. This 
discrepancy has led to some uneasiness that important information may be lost 
with V[ ow fc. Similar concerns have been expressed about effective field theory 
(EFT) potentials. 

This unease is exacerbated by experience with conventional NN potentials, 
which feature strong short-range repulsion. The repulsion causes even bound 
states with very low energies (such as the deuteron) to have important con- 
tributions to the binding and other properties from high-momentum com- 
ponents (well above 2fm~ ) of the deuteron wave function. In Ref. [Sj, the 
authors calculate cross sections for electron scattering from the deuteron us- 
ing as input a spectral function that is the momentum distribution times a 
delta function in energy. After noting the effect of excluding high momenta on 
the cross section, they conclude: "and thus the data confirm the existence of 
high-momentum components in the deuteron wave function" [8]. Beyond the 
fact that wave functions are not observables, these conclusions reinforce the 
intuition that there is information in quantitatively reproducing high-energy 
phase shifts that is lost when evolving to low-momentum interactions. In this 
paper, we use the similarity renormalization group (SRG) P|lUllllfl2] as a 
tool to demonstrate unequivocally that this intuition is incorrect. 

A fundamental tenet of renormalization theory is that the relevant details of 
high-energy physics for calculating low-energy observables can be captured in 
the scale-dependent coefficients of operators in a low-energy hamiltonian |13j . 
This principle does not mean that high-energy and low-energy physics is en- 
tirely decoupled in an effective theory. In fact, it implies that we can include 
as much irrelevant coupling to incorrect high-energy physics as we want by 
using a large cutoff, with no consequence to low-energy predictions (assuming 
we can calculate accurately). But this freedom also offers the possibility of 
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decoupling, which makes practical calculations more tractable by restricting 
the necessary degrees of freedom. This decoupling can be efficiently achieved 
by evolving nuclear interactions using RG transformations. (For an earlier dis- 
cussion of decoupling based on Okubo unitary transformations, see Ref. |14J.) 

The SRG allows for particularly transparent and convincing demonstrations 
of decoupling, because the evolution of the hamiltonian and other operators 
proceeds via transformations that can be chosen to be unitary, so that all 
observables are explicitly preserved. Thus, when we include the full set of mo- 
mentum states used for the original potential with the evolved potential V s , 
we find the same deuteron binding energy and phase shifts for all energies. 
However, the evolved V s explicitly decouples high-energy dynamics from low- 
energy observables, which means that we can exclude the high-momentum 
parts (so that we have a low- momentum potential like Mowfc) without disturb- 
ing the information content (phase shifts and deuteron binding energy) at 
lower energies. We apply a similar test to the SRG-evolved deuteron momen- 
tum distribution to show that high-momentum effects in low-energy bound 
states are captured by scale-dependent low-momentum operators. 



2 Background on the SRG 



In Ref. [12], the similarity renormalization group (SRG) approach is applied 
to NN interactions. The SRG was developed by Glazek and Wilson [9], while 
Wegner [10] independently developed a related but simpler set of flow equa- 
tions. The formalism we employ closely resembles that of Wegner, but we find 
that a transformation even simpler than advocated by Wegner is robust and 
adequate for all calculations to date. 

The evolution or flow of the hamiltonian with a parameter s follows from a 
unitary transformation, 

H s = U(s)HU\s)=T Iel + V s , (1) 

where T ve \ is the relative kinetic energy and H = T re \ + V is the initial hamil- 
tonian in the center-of-mass system. Eq. ([1]) defines the evolved potential V s , 
with T re i taken to be independent of s. Then H s evolves according to 

^ = [v(s),H s ], (2) 

with 

V(s) = ^U\s) = - V \s) . (3) 
Choosing r/(s) specifies the transformation. As in Ref. [12], we restrict ourselves 
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to the simple choice [TT] 

f](s) = [Trel, H s ] , 

which gives the flow equation, 



(4) 



= [[T rel , H S ],H S ] = [[T rel , • (5) 

In a momentum basis, this choice suppresses off-diagonal matrix elements, 
forcing the hamiltonian towards a band-diagonal form. 

The evolution in Eq. (pjj) includes all many-body components of the hamilto- 
nian. In the space of relative momentum NN states only, it means that the 
partial-wave momentum-space potential evolves as (with normalization so that 
1 = - Jo° < l 2 dq \q){q\ and in units where h = c = m = 1 with nucleon mass m), 



dV s (k, k) 2 r , 2 



ds 

2 r 00 



(k 2 - k' 2 ) 2 V a (k, k') 

O poo 

+-/ q 2 dq(k 2 + k' 2 -2q 2 )V s (k,q)V s (q,k'). (6) 

7T JO 



(The additional matrix structure of V s in coupled channels such as 3 Si- 3 D! 
is implicit.) For matrix elements far from the diagonal, the first term on the 
right side of Eq. ([6]) dominates and exponentially suppresses these elements as 
s increases. By discretizing the relative momentum space on a grid of gaussian 
integration points, we obtain a simple (but nonlinear) system of first-order 
coupled differential equations, with the boundary condition that V s (k,k') at 
the initial s = is equal to the initial potential. Since the SRG transformation 
is unitary, the NN phase shifts and deuteron binding energy calculated with 
H s are independent of s to within numerical precision. 

The evolution with s of any other operator O is given by the same unitary 
transformation, O s = U(s)OW(s), which means that O s evolves according to 

= fa( S )> O a ] = [[Trei, V s ),O s ] . (7) 

Just as with the hamiltonian H s , this evolution will induce many-body oper- 
ators even if the initial operator is purely two-body. If we restrict ourselves to 
the relative momentum NN space, we have 



dO s {k,k') 2 r°° 



poo 

/ q 2 dq (k 2 -q 2 )V s (k,q)O s (q,k') 
Jo 



ds 7T JO 

+ (k' 2 -q 2 )O s (k,q)V s (q,k>) 
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To evolve a particular O s simultaneously with V s , we simply include the dis- 
cretized version of Eq. ([8]) as additional coupled first-order differential equa- 
tions. 

An alternative, more direct, approach in the two-body space is to construct 
the unitary transformation at each s explicitly, and then use it to transform 
operators. Let \ip a (s)) be an eigenstate of H s with eigenvalue E a (which is 
independent of s). Then U(s) is given by 

u(s) = Y,\Ms))(Mo)\- (9) 

a 

In a discretized partial- wave relative-momentum space with momenta we 
can solve for the eigenvectors of if = H s= q and H s , then construct the matrix 
elements of U(s) [which we denote U s (ki, kj)] by summing over the product of 
momentum-space wave functions: 

u s (h,k J )=J2(h\Ms)}(Mty\k J }. (io) 

a 

In practice this is an efficient way to construct the unitary transformation and 
subsequently to evolve any operator in the two-body space. 



3 Applying the SRG 



In Ref. |12j . the SRG was applied with chiral EFT interactions at N 3 LO p31TT6] 
as the initial potentials. It was shown that the simple SRG transformation 
drives the hamiltonian towards the diagonal (in momentum space), making 
it more perturbative and more convergent in few-body calculations. (In fact, 
it behaves just like a Viowfc potential.) The parameter A = s -1 / 4 provides a 
measure of the spread of off-diagonal strength in V s . The V s potential at A was 
found to correspond roughly to a smoothly regulated Vi ow fc interaction [T7] 
with momentum cutoff A « A. In this work, we also use the Argonne fig [TS] 
and the CD-Bonn [B5] potentials as initial potentials in the SRG evolution. 
We find the same universal behavior [20] . 

In Fig. [TJ phase shifts in selected partial waves up to laboratory energies of 
1 GeV are shown for the Argonne vis [H], CD-Bonn [19] and chiral N 3 LO [15] 
potentials^ As expected for unitary transformations, the phase shifts at any 
A are the same as those for the corresponding initial potential (up to numerical 
inaccuracies that are currently less than 0.1% in the worst case, with no real 

1 Nonrelativistic kinematics is used, which means that comparisons with experi- 
ment at higher energy should be made with caution. This can be improved but will 
not change our discussion because it will have no effect on low-energy observables. 
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Fig. 1. Phase shifts for the Argonne vi$ [IB], CD-Bonn [19] and one of the chiral 
N 3 LO [15] potentials in selected channels (using nonrelativistic kinematics). The 
phase shifts after evolving in A from each initial potential agree for all A to within 
the widths of the lines at all energies. 
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Fig. 2. Phase shifts in selected channels for the Argonne vi% potential [18] and when 
intermediate momenta k > k max = 2.2 fm -1 are excluded. We contrast the latter 
results to the phase shifts obtained from the evolved V s potential for A = 2fm _1 
and the additional constraint k max = 2.2fm _1 . 
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attempt at numerical optimization). The phase shifts from different initial 
potentials start to disagree for energies above 300-400 MeV. The fact that 
they lead to very similar low-momentum interactions [TH27T2"] is an indication 
that all low-energy observables should similarly decouple after RG evolution. 
But here we explicitly demonstrate this decoupling rather than relying on 
generic RG arguments. 

We test the decoupling of high-energy details from low-energy phase shifts 
by setting V s (k, k') to zero for all k, k' above a specified momentum /c max . 
Results for /c max = 2.2fm _1 using a smooth regulator function are shown in 
Fig. [2] for the initial Argonne fig potential and for the evolved V s potential 
with A = 2.0 fm~ . The phase shifts for the initial potential in the lower 
partial waves bear no relation to the result without a £; max cutoff. The coupling 
between high and low momentum for the Argonne potential is so strong in 
the coupled 3 Si- 3 Dx channel that we had to use a slightly larger /c max = 
2.5 fm -1 with a smoother regulator just to keep the phase shifts on the plot. 
In contrast, the low-energy phase shifts for the SRG-evolved potential are 
unchanged, even though the high-energy phase shifts above fc max are now zero. 
The details of how the momentum is cut off affect only the behavior near /c max . 
The deviation near /c max in Fig. [2] is consistent with the SRG (with the present 
choice of rf) effectively imposing a rather smooth cutoff in momentum, similar 
to an exponential regulator with exponent n exp = 2, which is the analytic 
behavior for solutions of the linearized SRG equation. Finally, the angular 
momentum barrier in higher partial waves such as 3 F3 ensures that cutting off 
high momentum has the same effect for all interactions since they share the 
same long-range one-pion exchange potential. 

The deuteron binding energy provides another clear example of how the contri- 
butions of different momentum components to a low-energy observable depend 
on the resolution scale (as measured by A or s). In Fig. [3], we show the kinetic, 
potential, and total energy from an integration in momentum space including 
momenta up to fc max . That is, we plot 

Ed(k < A; max ) = 71 c i(A; < A; max ) + V s (k < fc max ) 

= J dkj dk'4(k; A) (k 2 5 3 (k - k') + V s (k, k')) A ) > 

(11) 

where ^(k; A) is the momentum-space deuteron wave function from the cor- 
responding potential V s (without k max ). 

Figure [3] shows that the Argonne vi$ potential [18] has large (and canceling) 
contributions from the high-momentum matrix elements of the hamiltonian. 
For example, if one excludes momenta greater than 2fm _1 in the deuteron 
wave function when calculating the binding energy, the deuteron is 9.9 MeV 
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Fig. 3. Expectation values in the deuteron of the kinetic, potential, and total energy 
evaluated in momentum space as a function of the maximum momentum fc max , see 
Eq. (fTTj) . Results are shown for the Argonne vis potential [18] (left), the evolved V s 
potential for A = 2fm~ 1 , and the smooth-cutoff V\ ow k interaction with A = 2fm _1 
and exponential regulator n exp = 2 

unbound (that is, the integrated kinetic energy up to 2 fm -1 is 11.5 MeV while 
the potential energy is —1.6 MeV). One needs to include contributions up to 
4fm _1 to even get a bound state. In contrast, using V s with A = 2fm~ , we 
see that the converged result is dominated by contributions from much lower 
momenta!" 2 ! Note that the V s potential has no appreciable contributions above 
A, even though the near-diagonal matrix elements of the potential V s (k, k') for 
k,k' > k max are not negligible. This again validates the advertised decoupling. 
We also see that the V s and V\ ow k results are very similar for A ~ A, where A 
is the momentum cutoff for V\ ow k- 

We turn to the momentum distribution in the deuteron next. The momentum 
distribution at relative momentum q is the expectation value of a^a^ (summed 
over spin substates M s ) and is proportional to the sum of the squares of the 
normalized S-wave and D-wave parts of the deuteron wave function, u(q) 2 + 
w(q) 2 . It is not directly related to an observable (see, for example, Ref. [21J). 
As discussed in Sect. [2[ using the SRG we can consistently evolve operators in 
s (or A). In particular, we can evolve a^a q starting from a given hamiltonian. 
Since the SRG proceeds via unitary transformations, no information is lost, by 

2 We can also simply set the high-momentum matrix elements of V s to zero (that 
is, exclude momenta larger than k max at the level of the hamiltonian), and then 
solve for the deuteron. For A = 2fm _1 , choosing & max = 2.2 fm" 1 only changes the 
binding energy to 2.2 MeV, while raising A; max slightly yields the full V s result. 
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q [fm" 1 ] 

Fig. 4. Deuteron momentum distribution (a q a q )d oc u{q) 2 + w{q) 2 using the Argonne 
Uis [15] . CD-Bonn [19] and SRG potentials evolved from Argonne vis to A = 1.5 fm -1 
and 2fm _1 (but not evolving the operator), and a smooth-cutoff Vjowfc interaction 
with A = 2fm~ 1 and exponential regulator n exp = 2. 

construction. But even more, we can show explicitly that by evolving to a low- 
momentum interaction, we decouple the low-momentum and high-momentum 
physics in a low-energy state. 

In Fig. HI we plot deuteron matrix elements of a q a q for the Argonne i>ig [IB] 
and CD-Bonn [T3] potentials, as well as for two SRG and a smooth-cutoff 
Viowfc interaction evolved from Argonne vi$. We emphasize again that matrix 
elements of this operator are not measurable, so one should not ask which of 
these is the "correct" momentum distribution in the deuteron; it is a potential- 
and scale-dependent quantity. It is evident that the V s distributions have sub- 
stantial momentum components only for k below A, and that the V\ ov ,k distri- 
bution is very similar to the corresponding V s distribution. Nevertheless, if we 
use the SRG- or RG-evolved operator with these deuteron wave functions, we 
precisely reproduce the momentum distribution from the original potential at 
all momenta and for all A. 

This result by itself is guaranteed by construction. The more interesting issue 
is where the strength of the matrix element comes from. For example, for the 
bare operator and the Argonne v is potential, the momentum distribution at 
q = 4fm _1 comes entirely from deuteron wave function components at that 
momentum. But when A = 2 fm -1 , it is clear from Fig.HJthat the deuteron does 
not have appreciable momentum components above 2.5 fm -1 (even though 
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Fig. 5. Ratio of the deuteron momentum distribution at various momenta q evolved 
from the Argonne v\g potential [18] via the SRG to the corresponding initial momen- 
tum distribution, as a function of the maximum momentum A; max in the deuteron 
wave functions in the numerator. Note that the un-evolved Argonne v\s result is 
simply a step function at q. For comparison, we also show the result for a smooth- 
cutoff V\ ow k interaction with A = 2fm~ 1 and exponential regulator n cxp = 2. 

V s does near the diagonal [12J). In Fig. El we take the ratio of the evolved 
operator evaluated with the evolved wave function at q to the corresponding 
initial quantity, but include in the numerator only momenta up to /c max . The 
numerator is thus: 



dk / dk'Vl(k;A)^(k,q)C7l(q,k')^(k';A) 



(12) 



We observe that for all q, the ratio approaches unity for large enough k max , 
as dictated by the unitary transformation. For A = 1.5 fm -1 or 2fm~ 1 , larger 
values of q (q > Sfm^ 1 or q > 4fm _1 , respectively) give results approximately 
independent of q, with a smooth approach to unity by & max m 1.3A. This is 
consistent with the operator U\(k, q) factorizing into K\(k) Q\(q) for k < A 
and q ^> A, and thus the q dependence cancels in the ratio. It remains to be 
seen whether this factorization is a general feature that can be understood 
using operator product expansion ideas [13]. For small q, the original step 
function is essentially preserved (not shown). For q of order A, there is a step 
behavior at A but some strength is shifted to lower k, k! . In all cases the flow of 
the operator strength weighted by ^ is toward lower momenta. The ratio for 
the smooth-cutoff V\ ow k is very similar at low fc max and for A; max approaching 
A/A. Finally, we have verified that taking matrix elements using even a rough 
variational ansatz for the wave function also works fine, which serves as a 
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check that there is no fine-tuning in the evolved operator. 



4 Conclusions 

Well-established renormalization theory tells us that low-energy physics does 
not depend on the details of the high-energy dynamics. The high-energy in- 
formation we do need can be incorporated in nuclear interactions using renor- 
malization group methods designed to handle similar problems in relativistic 
field theories and critical phenomena in condensed matter systems. This means 
that we can decouple the physics of low energy from high energy, drastically 
reducing the number of explicit degrees of freedom required for precise non- 
perturbative low-energy calculations. 

But while it is possible to decouple low- and high-energy physics, it does not 
happen automatically. We have seen that when using conventional NN po- 
tentials with strong short-range repulsion, such as Argonne vi% [H], there are 
momentum contributions that are high on nuclear scales but must be included, 
or else even very low-energy observables will be incorrect. We emphasize that 
the need for high-momentum components in these particular calculations does 
not imply that the high-energy description is correct (or measurable). 

Two-nucleon interactions derived via the SRG explicitly decouple low-energy 
physics from detailed high-momentum dynamics, allowing a clean demonstra- 
tion of these principles of renormalization. Using SRG potentials, we have 
shown that high-energy details in wave functions and operators are irrelevant 
to low-energy observables. When we set the high-momentum components to 
zero in V s potentials with low X, there is no noticeable loss of information: 
The low-energy phase shifts and expectation values in the deuteron are prac- 
tically unchanged. A corollary is that statements about high-momentum parts 
in the deuteron or any other bound state depend on the potential or scale. This 
also holds for expectation values of other operators, such as the one-pion ex- 
change potential or three-nucleon interactions. Finally, our results imply that 
there is also no problem with effective field theory interactions that predict 
zero phase shifts above a cutoff. 

Using the SRG, operators can be evolved, including those associated with high 
momenta in initial potentials, so that matrix elements in low-energy states 
(like the deuteron) are unchanged. We find for the momentum distribution 
that the operator strength flows toward lower momentum, so the matrix ele- 
ments do not change even after contributions from high momenta are excluded. 
We also observe that the ratio of evolved to initial deuteron momentum dis- 
tribution is nearly independent of momentum q at fixed A for sufficiently large 
q, which motivates an investigation using an operator product expansion. 
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In the past, unitary transformations of NN potentials were applied to cal- 
culations of few-body systems and nuclear matter. Results for observables 
depended on the transformations and this was often the context for discussing 
"off-shell" effects and which was the "true potential" . From the modern per- 
spective, this approach is misleading at best. These transformations always 
lead to many-body interactions, even if absent in the initial hamiltonian. While 
this fact was clearly recognized in past investigations [22], many-body forces 
were usually neglected. The SRG is one promising approach to evolve consis- 
tent three-nucleon interactions. Once the three-body SRG is implemented, the 
tests of decoupling can be extended to systems with more than two nucleons. 
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